data("ha.sts")
haagg <- aggregate(ha.sts, nfreq=13)
#A suite of of simple tests (inspired by the Matrix package)
stopifnot(identical(ha.sts, ha.sts[]))
plot(haagg[, 3])       # Single series
plot(haagg[1:30, 3])   # Somewhat shorter
#Counts at time 20
plot(haagg[20, ], type = observed ~1 |unit)Run the code above in your browser using DataLab